Decay and coherence of two-photon excited yellow ortho-excitons in Cu 2 0. 
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Photoluminescence excitation spectroscopy has revealed a novel, highly efficient two-photon exci- 
tation method to produce a cold, uniformly distributed high density excitonic gas in bulk cuprous 
oxide. A study of the time evolution of the density, temperature and chemical potential of the exci- 
ton gas shows that the so called quantum saturation effect that prevents Bose-Einstein condensation 
of the ortho-exciton gas originates from an unfavorable ratio between the cooling and recombination 
rates. Oscillations observed in the temporal decay of the ortho-excitonic luminescence intensity are 
discussed in terms of polaritonic beating. We present the semiclassical description of polaritonic 
oscillations in linear and non-linear optical processes. 
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10 ■ I. INTRODUCTION 

i— i' The quest for Bose-Einstein condensation (BEC) of excitons has been one of the main driving forces behind the 
numerous studies on high density excitonic gases during the last few decades^ A variety of systems have been tried, 
including bulk crystals^ 3 - and quantum well systems4^ In spite of the intense effort there still is no clear evidence that 
condensation was actually reached. The formation of a condensate is expected to occur if the Bose gas reaches a quasi- 
equilibrium state characterized by a sufficiently high density n and, simultaneously, a sufficiently low temperature 
T. For a three dimensional non-interacting Bose gas in equilibrium the transition temperature T c is proportional 
to n 2 / 3 . Strong spatial confinement and decreased dimensionality lead to a substantial decrease of the condensation 
temperature. This is why bulk crystals remain very attractive in the search for excitonic BEC. The long life times 
(> 1 ns) and high binding energy (150 meV) of the excitons in CU2O make cuprous oxide an excellent candidate in 
this class. In addition, the small Bohr radius (as ~ 7 A) of the ground state exciton assures that even at densities 
as high as n w 10 20 cm -3 the exciton-exciton interactions remain weak£ Apart from the intriguing possibility of 

11 BEC, there is a number of other interesting phenomena and properties of the exciton gas in CU2O which attracted 
attention, including localization of excitons in a stress trap^*£ spectroscopy on intra-excitonic transitions^ and the 

■ observation of polariton beating in transmission experiments resulting from interference of excited states belonging to 
different polariton dispersion branches^ii In this paper we revisit the high density ortho exciton gas in CU2O. The 
first part of the paper is concerned with spectroscopic photoluminescence excitation experiments revealing a novel 
and efficient non-linear exciton excitation scheme, which has several advantages over the schemes used so far. In 
ly-} , addition, we present a study of some of the statistical properties of the thus created exciton gas, and comment on the 
origin of the so called "quantum saturation effect"— that prevents BEC in the ortho-exciton gas in cuprous oxide. 
IO \ The second part of the paper concentrates on the polaritonic beating effect observed in our experiment for which 
a novel analytical model will be introduced. The new model explains both the beating reported earlieuifi^i for the 
transmission geometry as well as our results obtained for the back-reflection configuration. 

Cuprous oxide, CU2O, has four different excitonic series. These are attributed to various combinations of valence 
(V) band holes and conduction (C) band electrons (see Fig.^l: V\ — C\ (yellow), V2 — C\ (green), V\ — C2 (blue) and 
V2 — C2 (indigo) . The two uppermost valence bands originate from 3d copper orbitals (split by spin-orbit interaction) 
while the lowest conduction band is mainly formed from 4s copper orbitals. ■* Hence all three bands (Vi, V2 and C\) 
have the same even parity, and an optically forbidden direct gap (E s = 2.17 eV at 10 K) exists at a center of a Brillouin 
zone. The lowest excitonic state of the system is split by the electron-hole exchange interaction into a ls-ortho state 
(J = 1) at 2.033 eV (A ort ho) and a ls-para state (J = 0) at 2.021 eV (A para ). As a consequence of the even parity 
of the bands involved, the efficiency of "direct" optical pumping of the yellow excitons is very low. Radiative decay 
of Is ortho-excitons is weakly allowed through a quadrupolar coupling. For the ls-para-exciton radiative transitions 
are forbidden to high order. 

The yellow ortho-exciton emission spectrum consists of two main features: a direct k w transition (X™ ho ) at 
2.034 eV, and a phonon-assisted (A p r J ^ ho ) emission band at 2.02 eV (see FigflJ. The PA-band results from decay 
by simultaneous emission of a photon and a LO-phonon of energy 13.6 meV. Since both the matrix element and the 
phonon dispersion are nearly fc-independenlji^ the shape of the PA band directly reflects the kinetic energy distribution 
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FIG. 1: Left panel: Schematic energy band diagram at the F point in CU2O. Right panel: Ortho-exciton emission spectrum 
at T = 2.8 K using 532 nm continuous wave excitation. The thick line represents a fit to Eq.Q using a Maxwell-Boltzmann 
distribution with T bc = 24.7 K. 



of the ortho-exciton gas. The observed emission band may therefore be described by 

/oo 
D{x)f{x)e {S ^ L)2 dx , (1) 
-00 

where D{E) oc E 1 / 2 is the density of excitonic states (assuming a quadratic dispersion), and f(E) is either the 
Maxwell-Boltzmann (/ m b, classical gas) or Bose-Einstein (/be, quantum gas) distribution function. The product 
D{x)f(x) is convoluted with a Gaussian experimental resolution function of width T. Fitting this equation to the 
experimentally observed spectrum yields an estimate for the effective excitonic temperature and chemical potential 
(fi). For an ideal Bose gas above the condensation temperature, the density of bosons is given by 

n = g{m* /2nh 2 f' 2 J f he (E, n, T)D(E)dE , (2) 

with g the degeneracy, and m* the excitonic mass. 

The legitimacy of this analysis has been questioned in Ref. [lj, where it was argued that Auger processes might 
lead to a serious overestimation of the exciton density. Subsequent experimental studies^ however, have shown that 
this is only relevant near the surface. The experiments described here use excitation in the near infrared, where 
the optical penetration depth is larger than 500 (xmML Any contribution due to surface enhanced Auger processes is 
therefore expected to be negligible. For the same reason (large penetration depth) the thus created exciton gas in 
nearly uniformly distributed in the sample. Therefore effects related to a strong spatial inhomogeneity of the gas can 
also be neglected in our case. 



II. PHOTOLUMINESCENCE EXCITATION SPECTROSCOPY 



The samples used in this study were [100] platelets cut from a single crystal grown by a floating zone technique^ 
The upper limit for the full width at half maximum of the ortho-exciton direct emission at T=10K is measured to 
be 0.09nm (0.3mcV) which is an experimental resolution limited result. The excitonic emission is observed up to the 
room temperature. In the emission spectrum obtained for overbandgap excitation (E > 2.17 eV) a broad band is 
observed peaked at 1.7 eV attributed to ionized oxygen vacancy Vq + The relative intensity of the (A ort ho) and 
Vq + bands is comparable to those reported for high quality samples and shows that the samples used in this study 
have a relatively low density of the oxygen vacancies. Emission peaks attributed to excitons bound to impurities that 
have been reported around 2eV are absent in our spectral Moreover we observe an exceptionally long lifetime of 
paraexcitons in our samples, r w 14 ms, 19 independently demonstrating the high quality of the sample since the decay 
of paraexcitons is exclusively defect limited. 

For the photoluminescence excitation experiments samples were polished and mounted on the cold finger of a 
continuous flow cryostat and kept at T = 4 K. The time integrated experiments were performed using an optical 
parametric amplifier pumped by the third harmonic of a 10 Hz Nd:YAG laser. Excitation pulses with photon energies 
in the 1 — 2.5 eV range and temporal width of 8 ns were weakly focused on the sample (spot size 500 ^m). The 
yellow exciton emission was detected in a back-reflection geometry using a 0.5 m double-grating monochromator in 
combination with a photomultiplier. 



„ 3 

2 

C 

-° 1 



10 



1 wr 



1 10 
Power [mW] 



0.1 



detection window 



co 
b 



MB 



Mf, 



MB 
'../ 

oqd, 



42K 



JO 



R-TPE 



B-TPE 



OPE 



0.8 1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

£ [eV] 

pump L J 



FIG. 2: Integrated intensity of the ortho-exciton emission per incoming photon (IoE/n p ) versus the pump laser energy (i5 pump ). 
Bath temperature T — 7 K, P = 3 mW = constant. Full circles - Blue-exciton Two Photon Excitation (B-TPE), open circles 
- One Photon Excitation (OPE) and black dot - Is Resonant Two-Photon Excitation (R-TPE). Dashed lines are guides to the 
eye. Inset: Iqe as a function of the pumping laser power for -Ep Ump = 1.4 eV. Solid line shows a fit for I ~ P a with a — 1.8. 



Figure shows the results of a photoluminescence excitation study of yellow exciton creation in CU2O. During the 
experiments the pump fluency was kept constant at 30 W/cm 2 . The integrated intensity Joe of the ortho-exciton 
emission (both phonon-assisted and direct) is normalized by the number of incoming photons per second n p , and is 
therefore a measure of the efficiency of the pumping. 

Optical pumping is the most commonly used method to produce excitons. The efficiency of creating the excitonic 
gas will naturally depend on the particular process involved, and hence on the pump photon energy E pump . The 
observed excitation spectrum may be divided into several regions representing different excitation processes. A first 
rough division can be made into a one-photon (OPE) and a two-photon (TPE) excitation region, with the divider 
at the Is exciton energy, fn the high energy OPE region (i? pump > E gap w 2.17 eV) the dominant process is direct 
excitation of electron- hole (C± — V\) pairs, which may subsequently form an exciton. As mentioned above, however, 
this process is not very effective due to the even parity of both bands involved. For £ pump < E gap the process becomes 
more efficient since single photon transitions to excited yellow exciton states as well as resonant Raman transitions are 
allowed. The OPE efficiency reaches its maximum when E pump approaches A°j; ho -|-LO at 2.05 eV^i and decreases 
gradually upon increasing excitation energy. It agrees with the fact that the absorption coefficient at the excitonic 
states of higher quantum number n decreases as w 1/n 3 . Decrease of the number of free excitons created by pumping 
photon may also be related to the fact that overbandgap excitation enables retrapping of either electrons or hot- 
excitons by defect sitesi^l Moreover a surface photo-voltage effect may take place, which leads to a space-charge 
accumulation at the surface. 

Resonant two photon excitation (R-TPE, black data point at 1.015 eV) to the ortho-exciton state is allowed by 
symmetry and could thus, in principle, be efficient. Unfortunately, the matrix element for this process is fairly small 
since it involves an intermediate state which lies as deep as 5.6 eV below the upper valence band (Vi)jSi in agreement 
with the low efficiency observed. 

The so far unexplored part of the excitation spectrum is the 1.3-1.9 eV region. Surprisingly efficient excitation of 
the 2.036 eV yellow ortho-excitons is observed for E pump significantly below X ort ho but still well above A ort ho/2. The 
precise nature of the process responsible will be discussed later on. At this point, it is sufficient to note that it is a 
two photon process with Vi — C2 electron-hole pairs as intermediate state. This is evidenced by the observation of 
a threshold for the process at ~ 1.3 eV which coincides with half the V\ — C2 energy gap and by the near quadratic 
dependence of the excitonic emission on the fluency (see inset Fig.|2J measured for E pump =lAeV. 

The origin of the saturation of the efficiency for high fluencies (P » AQW/cm 2 ) remains to be investigated. It has 
also been observed for other excitation methods, and may originate from two-exciton spin-flip processes^ 

For the quest for Bose-Einstein condensation in excitonic matter it is vital to have an efficient excitation process 
(i.e. a high density). The limited lifetime of the exciton gas, however, also imposes the requirement of a low initial 
exciton gas temperature. At first sight, this seems problematic for the new B-TPE process. The process has an excess 
energy of at least 0.5 eV, which might be dissipated in either the exciton gas or in the lattice. A first estimate of 
the temperature of the gas can be made by fitting Eq.JTJ to the PA band (see Fig. QJ, using the Maxwell-Boltzmann 
distribution. The thus obtained gas temperatures are denoted in Fig. for several OPE and B-TPE excitation 
energies. In all cases the exciton temperature is higher than the bath temperature (4 K). In the OPE region, the 
smaller the energy mismatch between £ pump and A ort ho, the lower T m b. It decreases from 55 K for £ pump = 2.3 eV 



to 32 K for 2.15 eV. Surprisingly, an even lower gas temperature is observed when 1.4 < E puul p < 1.9 eV in spite of 
the fact that the energy mismatch (0.1-0.7 eV) is much larger than in the OPE cases. The exciton temperature in 
this range is found to be 29 K, nearly independent on the photon energy. Apparently, the excess energy leaves the 
system. One possibility would be that the V\ — Ci electron-hole state decays to the yellow exciton state through a 
dipole allowed transition either direct, or via an intermediate VI — CI electron-hole state. Another possibility is that 
the V\ — G2 state acts as an intermediate state for a hyper-Raman type process. In both cases the excess energy is 
emitted as a photon, which, in view of the low absorption coefficient at low energy, would indeed leave the sample. 
For determination of the precise nature of the process involved, one would need to perform emission experiments in 
the energy range 0.1-0.7 eV. 

III. TIME RESOLVED TWO PHOTON SPECTROSCOPY 

A more detailed insight into the properties of the B-TPE created exciton gas may be obtained from time-resolved 
luminescence experiments. Time resolved measurements were performed using a traveling wave optical parametric 
amplifier pumped by 1 kHz, 1.55 eV, 120 fs pulses produced by an amplified Ti-Sapphire laser system. The 150 fs 
excitation pulses (1.3-2.5 eV) were focused on the sample with a spot size of about 100/xm. The exciton emission was 
detected by a Hamamatsu streak camera system operating in photon counting mode (temporal resolution 10 ps). The 
experimental spectral resolution function is measured to be of a Gaussian type with r=0.37nm. 

Figure 13 shows the time-resolved ortho-exciton luminescence spectra after excitation with a 150 fs pulse (-E pU mp= 
1.4 eV) at T — 7 K. Initially, a broad spectrum is observed dominated by the phonon assisted transitions reflecting 
the relatively high excitonic temperature. Despite the rather indirect way of excitation, no initial growth of the 
luminescence has been observed, and the exciton creation appears to be instantaneous upon excitation within the 
experimental resolution (10 ps). As time progresses both the intensity and the shape of the emission change, where 
the former reflects the cooling of the exciton gas, and the latter the temporal decay of the exciton population. Direct 
and phonon-assisted peaks become well resolved for delay times longer than 0.2 ns, and show a similar long time 
decay dynamics. 
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FIG. 3: Time evolution of the ortho-exciton emission spectrum after excitation with a 150 fs, J5p Ump = 1.4 eV pulse. 

Two representative time-resolved A ort ho luminescence spectra are shown in figure 01 recorded 250 and 950 ps after 
excitation with a 150 fs pulse with E pnmp = 1.4 eV, together with fits to Eq.QJ, using a Bose-Einstein distribution 
function. The values for the gas temperatures at these times are found to be T = 34 K, and T = 17 K, respectively. 
As one can see, the fits reproduce the high energy tail very well. The deviations observed at the onset of the emission 
are due to replica-like structures which presumably appear due to a small strain induced by the mounting of the 
sample on the cold finger of the cryostati^ In the final fitting of the data this was taken into account by adding a 
scaled down replica of the spectrum, shifted by about 3 meV to Eq.QJ. 

The parameters resulting from fitting Eq. Q to time resolved spectra of the phonon assisted luminescence are shown 
in Fig. (upper panels). The two parameter sets shown correspond to experiments with excitation with 1.4 eV (filled 
symbols) and 2.17 eV (open circles), i.e. to B-TPE and OPE excitation, respectively. The cooling of the gas, shown 
in Fig. [Sji), shows a bi-exponential behavior. Initial fast cooling (r ~ 0.2 ns) occurs via optical phonon emission. 
Once the kinetic energy of the excitons becomes too low for this process, the relaxation occurs via acoustical phonon 
emission with a relatively long decay time (r ~ 6 ns). Though the cooling rate hardly depends on the excitation 
method, it is clear from the figure that B-TPE excitation yields lower gas temperatures due to the lower initial 
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FIG. 4: Left panel: Schematic energy band diagram at the T point in CU2O. Right panel: Time resolved ortho-exciton emission 
spectra at t — 250 ps and t = 940 ps. Dashed lines are a fit to Eq.Q using Bose- Einstein statistics with Tb c = 64 K; fi = 
-0.7 meV, and T bc = 23 K; fi = -0.24 meV, respectively. 

temperature. The chemical potential, shown in Fig. Eb), decays faster than the temperature, with a time constant 
t ~ 1.2 ns. This decay is essentially due to loss of particles. Consistently, the integrated intensity (Fig. Eh) of the 
luminescence shows a similar decay time (t ~ 1.4 ns). The conversion of ortho-excitons into para-excitons, the main 
cause of the loss of particles, strongly depends on temperature, hence the non-exponential behavior at early tiniest 
The observation that the particle decay is faster than the cooling rate, makes it implausible that the gas would undergo 
a BEC transition. This is demonstrated in the inset of Fig. Eh), which shows the (n,T) phase diagram. The state of 
the gas follows the BEC phase boundary, and no crossing is observed, similar to earlier observations^ Apparently, 
the gas adjusts its quantum properties to the density and temperature and is in quasi-equilibrium. 

As the gas cools down the relative occupation of the k w state is expected to increase. Indeed, the direct emission 
initially increases (see Fig. Eh), and the ratio of direct to phonon assisted emission increases with time (decreases with 
temperature, Fig. Eli)- For B-TPE excitation this ratio is higher than for OPE excitation for given Tbe- One can 
therefore conclude that the B-TPE process leads to a preferential occupation of the k ss state, which in itself is 
favorable for BEC condensation. 

Different explanations have been proposed for the so-called "quantum saturation" effect. In Ref. it has been 
suggested that the quantum saturation phenomenon originates from the fact that the ortho-gas is not in an equilibrium 
state so that the number of ortho-excitons is not well defined. As a consequence \x — by definition and no phase 
transition is possible. The present results do not support this conclusion. The kinetic energy distribution of the 
ortho-gas can be described with a very good accuracy by the BE distribution, and moreover decay of the chemical 
potential coincides with the particle decay obtained independently from the emission intensity decay. Therefore, one 
may conclude that the ortho-gas is in fact in quasi-equilibrium, and the particle density is a well defined number. 
Apparently the exciton gas adjusts its properties even though the phase boundary is not crossed. Kavoulakis 23 has 
argued that the above mentioned fitting based analysis to extract the exciton temperature and density is misleading 
since the ortho-gas is not in a quantum state but in a normal state and under this assumption the gas slides along 
adiabats given by n/riQ ~ 0.01, where uq is the quantum concentration. However, for this to be true, it would require 
an overestimation of the exciton density by two orders of magnitude, and it does not explain how a classical gas can 
show a non-classical kinetic-energy distribution. A few possible explanations for an apparent non-classical distribution 
have, however, been given in literature: dominant Auger recombination and/or strong spatial inhomogeneity. As has 
been discussed in section I, both effects are expected to have only a marginal influence on the decay kinetics for B-TPE 
excitation. Since for OPE pumping similar results are found, the significance of these two effects is questionable in 
that case as well. Finally, Ell et alM presented an analysis of the ortho-exciton relaxation which included polaritonic 
effects. The main conclusion of that analysis is that the so called polaritonic "bottleneck effect" prevents formation 
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FIG. 5: Upper panels: Time dependence of (a) the exciton temperature Tb e and (b) chemical potential. Solid lines are fits to a 
two- and one-exponential decay, respectively. The inset in (a) shows the n, T phase diagram (solid line is the BEC boundary), 
where n has been calculated using Eq.(j2j. 

Lower panels: (c) Temporal decay of the integrated intensities of the total, direct (DI) , and phonon assisted (PA) luminescence. 
Solid line: fit to a single exponential decay, (d) Ratio between direct and phonon assisted luminescence as function of the gas 
temperature, solid lines are a guide to the eye. The data sets are for -Ep Ump = 1.4 eV (circles), and E puTnp = 2.17 eV (squares). 



of an ortho-excitonic BEC at k=0. However, as pointed out in Ref. l2?l the authors assumed in their model an infinite 
volume and it is not clear to what extend their conclusions would change due to finite volume corrections. From the 
above analysis we can conclude that only the third model does not contradict with the present results but the full 
model is still to be established. 

Aware of the limitations mentioned in the introduction, we applied the same analysis to £ , pump =2.17 eV data, 
where the penetration depth is only 30 /im, and the results lead to similar conclusions as in the B-TPE case. The 
fact that the early time density reached for £ , pump =2.17 eV excitation is significantly higher than for _E pump =1.4 eV 
excitation agrees well with the ~10 times smaller penetration depth, and the 5 times lower pumping efficiency. This 
indirectly supports the use of Eq.©, even in the case of a small penetration depth. 



IV. POLARITON OSCILLATIONS 



Analysis of the X ort ho emission presented in the previous sections shows that absorption of the pump photons 
results in the formation of ortho-exciton through a fast and efficient process. The intensities of both X^ ho and 
^ortho emission decay mainly as a result of down-conversion of the ortho-excitons into para-excitons. More detailed 
analysis of the X™ ho however reveals also weak oscillations in the measured X^ ho emission intensity (I d), reminiscent 
of those previously observed by Frohlich ct al. 10-11 This is shown in Fig. which displays the time dependence of the 
intensity after subtraction of a smooth exponential decay curve. The oscillations have a typical "period" of about 0.5 
ns, which increases with time as y/t (see also inset figEJ- They are thought to originate from interference between 
states of the lower and upper polaritonic branches i&iL Such oscillations have so far only been observed for resonant 
one-photon excitation in a transmission geometry, and it is not immediately clear that the polaritonic interference 
would be relevant to the present two-photon experiment and a B-TPE excitation scheme. In the following section we 
therefore make an excursion to the semiclassical theory of the polariton oscillations of transmitted light, which allows 
us to obtain an analytical expression for these oscillations. Most importantly, the results obtained in Sec. |V] can be 
generalized to describe more complicated processes, such as the present time resolved two photon experiment. This 
is done in Sec. IVII where we show that the aperiodic oscillations presented in Fig. H3 indeed result from the polariton 
beating. 
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FIG. 6: Direct emission versus time from figure 5(c) after subtraction of a single exponential decay curve (symbols). Dashed 
and solid lines: theoretical simulation as discussed in sectior TVTl The inset shows the \/t dependence of the peak positions. 



V. POLARITON OSCILLATIONS IN TRANSMISSION OF LIGHT AND GEOMETRICAL OPTICS 



In this section we give a semiclassical description of the polariton oscillations, observed in the transmission experi- 
ment on Cu^O^ii In principle, the transmission intensity can be found by a straightforward numerical calculation, 
as was done in Ref. flfT ] . Here we obtain an accurate analytical expression for this intensity. 

Consider the propagation of linearly polarized light with frequency u through a slab of a material with dielectric 
function e (see Fig. [7J) . For a monochromatic incident electric field Ei = Eoe~ luJt+lkx , where k = — , the amplitude of 
the transmitted field is 

E t= E oe- iUt+iq \^&% > (3) 

where L is the slab thickness, q = ^\fe, and R = ^fjj^; is the reflection amplitude. Equation @ allows for a simple 
interpretation: the electric field is multiplied by the factor 1 — R when the light enters the material, and by 1 + R 
when it gets out, — itot + iqL is the phase acquired by the light traveling through the sample during the time t, and 
the denominator 



1 - R 2 e 2l i L 



1 + i?-V^ + We 



is the sum over the number of internal reflections from the boundaries the slab (see Fig. |HJ: a factor R for each 
internal reflection and e %qL for the propagation from one side of the slab to the other side. 
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FIG. 7: The transmission of light propagating in the x direction through the slab. 



Consider now the incident electric field having the form of the Gaussian wave packet with the average frequency u 
and width A: 

E in = / due 2^ . 

V2^A J 
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n = l 



71=2 



FIG. 8: Examples of processes with different numbers of internal reflections. 



Then the transmitted field at x = L is given by 



E * = -#jr E t d " i 1 - R2 ) ^ 2 v s "-^#, (4) 



n— 



where the action S n = —u>t + q (2n + 1) L is the phase of the field for a path with 2n internal reflections. Since 
L is usually much larger than the wavelength, the integral in Eq.Q can be calculated in the stationary phase 
approximation: 



fl - R 2 ) R 2n 



E t ^E ^ ^ "l" e ^-^, (5) 

„=o Jl-i(2n + l)LA^ 



where the values of S n , -j-f , and R are taken at the 'stationary' frequency w n , such that 



-K)+^ = -(|-(2n|l)L + lT = 0. (6) 



For a broad wave packet, the term oc can be neglected and we obtain 

du _ (2n+ l)L 

dq 



(7) 



which is the geometrical optics result: the group velocity of light in the slab equals the path (2n + 1) L divided by 
the time of motion t. Furthermore, for a broad wave packet 



where A„ = 



E t « E £ (1 - i? 2 ) iZ*v( s »+* ), (8) 

-1/2 



A 

n=0 



(2n + 1) L4-i <C A is the width of the frequency interval around oj n that gives the dominant 

contribution to the integral in Eq. ( QJ . 

If there are no excitations with u) ~ u> present in the sample, e is a slow function of frequency and can be replaced 
by a constant eo = e o + * e o- The time it takes light to traverse the sample is then to = 7\/eo (f° r a l mm slab of 
CU2O with e = 6.5, to — 8.5ps). If the time t between the moment the light enters the sample and the moment it 
is observed on the other side of the sample greatly exceeds to, then one way to satisfy Eq.(JZJ) would be to consider 
processes, in which the light is reflected many times (n ^> 1) and passes a long way (2n + 1) L before it leaves the 
sample. Due to the rather strong absorption of visible light in CU2O and the loss of intensity at each reflection (the 
reflection coefficient R 2 w 0.2), the intensity of the transmitted light for n ^> 1 is extremely weak. 

Equation Q can also be satisfied for small n, but only if light in the sample can propagate slowly. This happens 
when it is mixed with an exciton, i.e. when Q ss f2, where fl is the exciton frequency. Such a mixing results in 
two polariton branches w± (q), both containing slowly propagating excitations (see Fig. EJ, so that Eq.Q has two 
solutions. Then Eq.© has to be modified and for each n the transmitted field is the sum of the contributions of the 
two polariton trajectories: 

n=0 ^ ' 



q 



FIG. 9: Two branches of the polariton dispersion (solid lines). The circles mark a pair of polaritons moving with the same 
velocity. 



where ± labels the polariton branch. Below we show that for the quadrupolar of excitons to light 5 " + 2 5 " = 
(2n + 1) ftto, where / is the strength of the quadrupolar coupling, so that the interference between the two 
polaritons moving with the same velocity gives rise to oscillations of the intensity of the transmitted light with a 
time-dependent period. Equation @ is actually an asymptotic expression for E t valid in the large-i limit, when 
Lu n+ — u n - A„, i.e. when the intervals of relevant frequencies around the two stationary points do not overlap. A 
more general expression is obtained below [see Ea. (|16|l ]. A similar situation is found in scattering problems, when the 
interference between different classical trajectories makes the scattering cross-section a rapidly oscillating function of 
angle. A simple semiclassical expression for the cross-section is not valid, however, for angles close to and tt, where 
the so-called glory phenomenon takes placed 
The polariton dispersion is found from 

(uj 2 + iloT - e 2 ) (uj 2 - uj 2 q ) - fcoW q = 0, (10) 

where uj q — and e q = Q + j^- are, respectively, the photon and exciton frequencies, / is the strength of the 
quadrupolar photon-exciton coupling and T is the exciton decay width. Equation (|10|) corresponds to the dielectric 
functioniSiii 

fc 2 q 2 

e(uj,q) = e + — - — . (11) 

e q — uj z — iuil 

2 2 2 

It is convenient to introduce the dimensionless variables X — f^r (wave vector squared) and Y = (frequency 

squared), which near the crossing point of the photon and exciton dispersions are both close to 1. We also make use 
of the fact that for CU2O the dimensionless parameters, describing the exciton-photon coupling (/ ps 4- 10 -9 ), exciton 

kinetic energy (£ = e' ^ sa 1CT 5 ), decay width (7 = £ sa 5 • Vd~ T )?ML and the absorption of light (c = ^ ~ 10~ 4 

are rather small. Then Ea. (|10fl can be written in the form 

(Y-l-Z + i 7 )(Y-X+i() = f. (12) 

Near the photon-exciton crossing point 

X = 1+Z + X', 

y = i + e + y, 

with X',Y' 1. Using Eq. l|12|l we can express the polariton 'momentum' X' in terms of its 'frequency' Y': 

X^Y'-yJ^+iC (13) 

Note that while X' (Y 1 ) is a single-valued function, Y' (X 1 ) has two branches. 
The polariton action is then given by 

S n = -ujt + q (2n + [(2n + 1) *o - *] + y [r n X' -Y'-£], (14) 



where r n — [2n + 1) % and, by assumption, r n <C 1. Using the substitution Y' = z^/fr n and Ea. (|13[l . we obtain 



a w c-(O) _ 9n 



z + 



z + i- 



(15) 



Here S*i 0) = fi [(2n + 1) t - t] + in (2n + 1) i J is the action for / = and <j> n = Qty/fr^ = Qy/(2n + 1) ftt is the 
phase of the polariton oscillations. Equation (Jlj can now be written in the form 



E t KE (l-R 2 )Y,^R 2n e iS ^l(<l>n 



n=0 



7 <t>n i 



(16) 



where Al = ^4^, i? = ^3U, and 

V4+ 1 



4t 



/(a, 6, c, d) = 



+ OC 



(17) 



The factor e 2( z + d ) (with c = (^f^) 2 and d = £ ), describing the shape of the Gaussian wave packet with a; = O, 

is only important for a < 1. For a» 1, I(a,b, c, d) becomes independent of c and d. In this limit the integration 
over z can be done in the stationary-phase approximation, as the two stationary points z n ±: z n ± + i/3 n = ±1, are 
sufficiently separated from each other. The result is 



I(a, b, c, d) ~ 2e ^ cos (a + ^ 



(18) 



In Fig.^]we compare I (a, b, c, d), defined by Ea. H17|) . for b = c = 0.1 and d = (dots), with the asymptotic expression 
Ea. (|18J) (solid line). We see that after about one oscillation the two curves coincide and the transmitted field equals 
the sum of the two polariton contributions 



E t « 2E (1 - R 2 ) ]T fl2n e in[(2„+i) t0 -i]-i Qn (2„+i)i-J? cos fn- s /(2n + l)ftt + - 



71=0 



(19) 



[cf. Eq.@], where a(fl) — -^7= is the absorption coefficient at the frequency O. 




FIG. 10: The a-dependence of the integral /(a, 0.1, 0.1, 0) [see Ea. l|17|l ] calculated numerically (dots) and in the stationary 
phase approximation (solid line). 



For further discussion it is useful to note that Eq. i|19fl can be written in the form 



(20) 



n=0 



where F (l,t) is a function of the length of the polariton path I and the time of motion t: 

■i«(n)i- 

F(l,t) = 2e 
l n = (2n+l)L and A (I) = A 



(21) 




FIG. 11: The ratio of the intensities of the transmitted and incident light for the parameters taken from Ref. [TflfTlj ] calculated 
numerically using Eq. H6H (dots) and in the stationary phase approximation (solid line). 



The ratio of the intensities of the transmitted and incident light for the parameters taken from Ref. [ToBTl| is shown 
in Fig. ^2 The results of the numerical calculation of Ecis. l|16l) and {T7J|] are plotted with dots, while the solid line 
was obtained using the semiclassical expression Eq. (|19|) ■ The difference between the two is hardly visible in this semi- 
logarithmic plot, even for short times ~ 30 ps. For the rather thick sample used in that experiment [L = 0.91mm) 



the contributions of the processes with multiple internal reflections are negligibly small: 
the transmitted intensity shows a single-mode oscillation, the period of which grows with time. 



h(n=l) 
h(n=0) 



2 • 10" 3 , so that 



VI. POLARITON OSCILLATIONS IN NON-LINEAR PROCESS 



In this section we return to the present experiments and apply the semiclassical formalism discussed in the previous 
section to a non-linear process, in which two photons with the frequency luq are converted into a photon with the 
frequency close to the ortho-exciton frequency Q < 2loq. Our aim is to explain the oscillations of the emitted light 
shown FigEl Although small in amplitude and somewhat irregular, these oscillations occur on the same time scale 
(of a few hundred picoseconds) as the oscillations of the transmitted light, discussed in Sec0 Furthermore, the inset 
in FigH3 shows the phase of these oscillations grows proportionally to square root of time, just as in the transmission 
experiment. Clearly, the intensity oscillations observed in our experiment result from the polariton beating. 

This conclusion is actually very surprising, since the polariton effects seem to be incompatible with the high-density 
exciton gas, achieved in the present set-up. The exciton collision time in a gas with the density n = 10 19 cm -3 and 
temperature T — 30K is r ~ 0.06ps, which is three orders of magnitude smaller than the oscillation period ~ 500ps. 
Furthermore, the polariton momentum p p — y^^r is much smaller than the width of the thermal distribution in the 
exciton gas pr = ySMT: p p ~ O.lpT- Thus the elastic collisions effectively suppress the polariton oscillations by 
removing very fast the excitons from the momentum range, in which their interaction with light is important. The 
intensity oscillations may only come from polaritons propagating through the regions where the exciton gas density 
is sufficiently low. In our experimental set-up the incoming light falls at a small angle ~ 10°? with respect to the 
surface normal, while the aperture angle was made rather large ~ 20°? to increase the signal-to-noise ratio. This 
allows polaritons, created at the rims of the light spot, to reach detector without passing through the dense exciton 
gas regions. The fact that the exciton cloud grows with time does not alter this conclusion. While the polariton 
momentum is much smaller than the thermal momentum, the polariton velocity v p ~ ~^z^t' where t ~ 2ns is the 

time of the polariton motion, is much larger than the thermal velocity vt = \J ^j-: v p ~ 20vt- Furthermore, the 
length of the polariton path ~ L covered during that time is much larger than increase of the size of the exciton cloud 
~ y 6vj.tr, so that polaritons can escape the expanding exciton cloud. 

The essential difference between the non-linear and transmission experiments is the energy and momentum dissi- 
pation, which occurs when the two in-coming photons are converted into a slowly propagating polariton. As details 



of this process are poorly understood, we simply assume here that the conversion is fast (on the scale of the polariton 
oscillations) and that it occurs locally, i.e. the polaritons are almost immediately emitted from the spot, where the 
two photons were absorbed. Furthermore, we assume that after the short thermalization time the motion of slow 
polaritons (with the frequency close to £1) is essentially free, i.e. their interactions with phonons and other excitons 
can be neglected. In the transmission experimentpiSiii the frequency of the incident light was deliberately chosen to be 
close to the exciton frequency f2 to maximize the amplitude of the polariton oscillations. In the nonlinear experiment 
one obtains an exciton gas with a broad thermal distribution. The energy and momentum selection of the excitons 
occurs in this case automatically, since only the excitons with the energy close to hfl are strongly mixed with light and 
contribute to the polariton oscillations, and only those moving in the directions approximately parallel and antiparallel 
to the sample surface normal have a chance to be detected (examples of such processes are shown in Fig ll2|l . This 
selection, however, strongly reduces the relative magnitude of the polariton oscillations. A further reduction comes 
from the fact that the energy and momentum dissipation destroys the interference between the polaritons created at 
different points in the sample. In other words, a luminescent exciton interferes only 'with itself. It is surprising that 
even in this case the quantum polariton oscillations are not entirely suppressed and can be detected. 




£1 < 

L L 



W (d) 



FIG. 12: Examples of nonlinear processes, in which two photons with frequency coo are converted into a polariton with 
frequency Q. For clarity, the directions of propagation of the incident and emitted light are drawn slightly non-parallel. 

In our simple model the intensity of the emitted light can be calculated as follows. Using the semiclassical method 
discussed above, we first find the radiation from the polaritons created at a point with the coordinate x in the sample. 
Then this intensity is multiplied by the conversion probability, weighted with a factor describing the absorption of 
the incident light in the sample, and integrated over x from to L. 

Consider, e.g., the process shown in Fig. 112b . in which two incident photons with frequency ujq enter the sample 
(with the probability (1 — i?) 4 ), the polariton travels along the path of length x in the direction opposite to that 
of the incident light, which gives the factor F 2 (x,t) [cf. Ea. (|20fl ] . and then leaves the sample (the corresponding 
probability is (1 + -R) ). For simplicity we assume the dielectric function for / — to be frequency independent, so 
that the reflection of both incident and emitted photons is described by the same coefficient R. The intensity of this 
process is then given by 

L 

J« oc (1 - Rf (1 + Rf J dxe~ 2a{uJa)x F 2 (x, t) . (22) 
o 

Here the factor e - 2a ( u o)x describes the suppression of the intensity of the incident light squared (a (wo) « ^a(fl)). 
Actually, the time of the polariton motion is somewhat shorter than the total time t between the moments light enters 
and leaves the sample. However, the conversion is, by assumption, fast and the time of motion of the incident photons 
in the sample of thickness L = 0.3mm, used in this experiment, is only a few picoseconds, whereas the oscillation 
period is of the order of one nanosecond. 

In the process shown in Fig. 112b . the polariton continues to move in the direction of the incident light and is 
reflected from the right boundary before leaving the sample (the path length is 2L — x) . The reflection results in the 



extra factor R 2 : 



L 

I b oc (1 - i?) 4 R 2 (1 + Rf J dxe- 2aMx F 2 (2L - x,t) . (23) 

o 

Similarly, the sum of the intensities of the processes shown in Fig. 11 2h and d is 

L 

h + I d oc (1 - i?) 4 i? 4 (1 + R) 2 J dxe- 2aM{2L - x) [F 2 (x, t) + R 2 F 2 (2L - x, t)] . (24) 

o 

In the same way one can find the emission intensity for the processes with multiple internal reflections of the incident 
and emitted photons. 

This simple approach gives a reasonable description of the observed intensity oscillations, as shown in Fig.0 where 
we plot the oscillating part of the intensity, I osc — Ie rt , measured experimentally (circles) and calculated theoretically. 
The dashed line is the sum of the contributions of all four types of processes depicted in Fig.^|with multiple reflections 
included. We noted that the agreement with experiments improves, if the contributions of the processes shown in 
Fig. 1121 a and d, in which the polariton is emitted in the direction opposite to that of incident light, are suppressed 
by the factor of 3. The rationale for doing this may be the fact that not all polaritons can reach the detector, since 
some of them are absorbed in the regions with dense exciton gas. Therefore, each of the processes discussed above 
has to be multiplied by a geometrical factor equal the fraction of excitons that can be detected. We did not perform 
the extensive calculations of this factor, since the results shown in Fig. El clearly show that the intensity oscillations 
result from the polariton beatings. 



VII. CONCLUSION 



In conclusion, the here presented B-TPE process provides an efficient method to create a high density, cold ex- 
citon gas in the bulk of cuprous oxide. The low initial temperature of the gas (despite the energy mismatch), the 
instantaneous observation of X ort ho luminescence, and the preferential occupation of the k sa state seem to point 
to a coherent type of B-TPE process as opposed to simple two photon absorption, and subsequent incoherent decay 
to the ortho exciton state. One possible process which would at least explain the low initial gas temperature would 
be two photon absorption from the upper valence band (^7) that excites a continuum of excitonic p states associ- 
ated with the upper conduction band (~T S )*& Due to the spin selection rules the excitons created this way will be 
triplet excitons. However the C2 band originates from p orbitals so we do not have pure spin states for the electron 
anymore. Therefore the angular momentum selection rules allow the Ci triplet exciton to decay radiatively (dipolc 
allowed transition) to cither triplet or singlet C± exciton and this way both para and ortho states of a yellow series 
may be populated. This could actually even happen in a coherent process (hyper- Raman scattering), with the V1-C2 
electron-hole pair as the intermediate state. The details of the B-TPE process remain unclear at the present. Future 
experiments on photon emission coming from the C2-C1 transition should yield the necessary information. 

We gave simple analytical expressions for the polariton oscillations observed both in the transmitted light in case 
of the resonant excitation and in the emitted light in time resolved two photon experiment. While in the first case 
the result is obtained by a straightforward application of the semiclassical method, the description of the beating that 
originates from a non-linear process requires additional assumptions. Since we have not found any evidence for the 
coherent production of the ortho-excitons, we assumed that the coherence of the incident light is lost in the irreversible 
conversion of two photons into a low-energy exciton. We argue that this does not completely suppress the polariton 
oscillations, since some polaritons can avoid the areas of dense exciton gas on their way to detector and live long 
enough to be considered as free particles. After such polaritons are created, they propagate as quantum superpositions 
of two polaritons, which gives rise to quantum beating, similar to the strangeness oscillations in K-mesons produced 
in hadron reactions .[A. Angelopoulos et al., Phys. Lett. B 503, 49 (2001).] The shape of the oscillating intensity of 
light obtained within this simple approach is in good agreement with the experimental data. 

Although the present study shows that ortho-exciton BEC can not be achieved in this kind of experiments, mainly 
because the particle decay rate is faster than the temperature decay, the TPE method does provide an efficient method 
to create a cold exciton gas. In particular, it allows to create a para-exciton gas in conditions where ortho-para down- 
conversion is blocked (i.e. at low temperatures), opening new possibilities in the excitonic BEC quest. 
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